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PREFACE 

This  report  covers  work  carried  out  as  a  joint  program  between 
engineers  from  the  Aero  Propulsion  Laboratory  (AFWAL/POTP)  and  the 
National  Aeronautics  and  Space  Administration  (NASA)  Lewis  Research 
Center  under  in-house  project  30661252.  The  objective  of  this  effort 
was  to  analytically  determine  the  elastic  stress/strain-temperature-time 
history  at  the  critical  location  for  a  double  edge  wedge  geometry 
specimen  cycled  in  fluidized  beds. 

The  research  was  conducted  from  June  1977  to  January  1979. 
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SECTION  I 
INTRODUCTION 

One  important  area  of  research  necessary  for  advancing  the  technology 
of  aircraft  gas  turbine  engines  is  the  accurate  assessment  of  the  life 
prediction  procedures  used  for  hot  section  blades  and  vanes.  In  order  to 
further  develop  and  evaluate  life  prediction  methods,  this  program 
tested  in  the  laboratory  simulated  hardware  components  using  carefully 
controlled  conditions.  Comparison  of  the  experimentally  measured  life 
to  that  which  is  analytically  predicted  Is  used  as  a  means  of  evaluating 
life  prediction  theories. 

The  experimental  laboratory  method  used  in  this  program  for 
measuring  thermal  fatigue  life  is  the  cycling  of  wedge  (blade-like) 
specimens  in  fluidized  beds.  Such  tests  have  been  shown  to  provide  life 
and  transient  temperature  data  under  carefully  controlled  conditions. 
Reference  1  contains  a  compilation  of  such  data  including  a  description 
of  the  facility  and  test  procedure.  References  2-9  contain  incremental 
portions  of  such  data  relative  to  the  evaluations  described  in  this  paper. 

The  objective  of  this  investigation  was  to  analytically  determine 
the  elastic  stress/strain-temperature-time  history  at  the  critical  location 
for  a  double-edge  wedge  geometry  specimen  cycled  in  fluidized  beds.  This 
was  performed  as  a  joint  program  between  the  engineers  from  the  Air  Force 
and  the  National  Aeronautics  and  Space  Administration  (NASA)  and  utilized 
conventional  three-dimensional  finite  element  elastic  analysis  techniques. 
Engineers  at  the  Aero  Propulsion  Laboratory  (AFWAL/POTP)  used  the  IS03DQ 
computer  program  while  the  NASA/Lewis  Research  Center  engineers  used  the 
NASTRAN  program.  The  alloys  were  IN  100,  Mar-M  200,  War-M  302,  NASA 
TAZ-8A,  and  Rene  80. 

Two  fluidized  beds  were  used  for  rapidly  heating  and  cooling  the 
specimens.  The  specimens  were  in  the  form  of  prismatic  bars  with  a 
double-wedge  constant  cross-sectional  geometry.  These  specimens  failed 
by  thermal  fatigue  cracking  which  is  usually  the  predominant  failure 
mode  of  aircraft  engine  first  stage  turb1ne_blades  and  vanes.  Thermal 
fatigue  is  defined  as  the  cracking  of  a  material  Induced  from  cyclic 
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stresses  and  strains  caused  by  repeated  temperature  changes.  The  cycling 
condition  was  alternate  3  minute  immersions  in  fluidized  beds  maintained  at 
316°  and  1088°C  (600°  and  1990°F).  The  cycling  test  condition  chosen  was 
one  which  resulted  in  thermal  fatigue  cracking  in  a  reasonable  number  of 
cycles. 

Due  to  symmetry,  a  discretized  model  of  only  a  quarter  of  the  double¬ 
edge  wedge  geometry  was  necessary  for  analysis.  First,  a  model  with  a 
fine  mesh  for  IN  100  alloy  for  a  severe  time  increment  (15  seconds  after 
immersion  in  the  heating  bed)  was  analyzed  using  the  NASTRAN  computer 
program.  Then,  a  model  with  various  coarse  meshes  for  the  same  conditions 
was  analyzed  using  the  IS03DQ  program.  A  coarse  mesh  model  for  the 
IS03DQ  analysis  was  selected  which  gave  essentially  the  same  results  as 
using  the  fine  mesh  model  with  the  NASTRAN  analysis.  The  remaining 
combinations  were  then  analyzed  using  the  coarse  mesh  model  and  the 
IS03DQ  program.  Such  analyses  provide  the  strain  range  and  stress/strain- 
temperature-time  history  so  important  for  evaluation  of  life  prediction 
theories.  The  turbine  component  life  prediction  methods  currently  being 
studied  at  NASA/Lewis  are  discussed  in  References  10  -  15.  Results  of 
similar  analyses  for  a  single-edge  wedge  geometry  specimen  are  given  in 
Reference  16. 
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SECTION  II 

INPUT  FOR  COMPUTER  PROGRAMS 

The  alloys  and  test  condition  for  the  five  alloys  are  given  in 
Table  1.  The  necessary  inputs  to  perform  the  analyses  were:  (1)  the 
geometry  of  the  double-edge  wedge,  (2)  the  elastic  and  physical  material 
properties  of  the  five  alloys,  and  (3)  a  complete  temperature  distribution 
at  various  times  throughout  the  cycle.  This  section  gives  a  detailed 
description  of  these  inputs. 

1.  WEDGE  GEOMETRY 

The  geometry  for  the  double-edge  wedge  is  shown  in  Figure  1(a).  The 
computer  plots  of  the  models  used  for  analysis  and  a  typical  element  are 
shown  in  Figure  1(b)  for  the  IS03DQ  program  and  in  Figure  1(c)  for  the 
NASTRAN  program.  The  model  for  both  programs  "squared-off"  the  leading 
edge  radius  to  a  1.02  mm  (0.040  in.)  length  and  the  trailing  edge  radius 
to  a  1.53  mm  (0.060  in.)  length.  Otherwise  the  models  duplicated  the 
geometry  of  the  wedge  exactly.  Detailed  discussion  of  the  modeling  is 
given  in  Section  III  Description  of  Analyses. 

2.  ALLOY  PROPERTIES 

The  temperature  independent  and  temperature  dependent  alloy  properties 
used  for  the  elastic  analyses  are  given  in  Tables  2  and  3,  respectively. 

The  properties  required  for  the  analyses  were  Poisson's  ratio,  modulus  of 
elasticity,  and  the  mean  coefficient  of  thermal  expansion.  The  programs 
required  a  value  for  density  to  obtain  results  (zero  mass  elements  were 
not  permitted)  although  the  results  are  independent  of  density.  The 
properties  for  all  alloys  except  the  mean  coefficient  of  thermal  expansion 
for  NASA  TAZ-8A  alloy  were  obtained  from  References  17  and  18.  The  mean 
coefficient  of  thermal  expansion  for  NASA  TAZ-8A  was  independently 
determined.  This  and  all  data  in  Reference  17  were  determined  from  the 
same  heat  used  for  fabricating  the  double-edge  wedge  test  and  calibration 
specimens. 


3 


AFWAL-TR-80-201 3 


3.  TEMPERATURE  LOADING 

The  transient  temperature  loading  on  the  double-edge  wedges  was 
determined  from  thermocouple  data.  Calibration  specimens  of  the  five 
alloys  were  instrumented  chordwise  at  the  mid-span  with  five  embedded 
thermocouples  and  cycled  in  the  fluidized  beds  (schematically  shown  in 
Figure  2).  The  location  of  the  thermocouples  at  the  wedge  cross-section 
is  shown  in  Figure  3.  The  Inconel  600  sheathed  thermocouples  were  mounted 
in  grooves  milled  in  the  surface  of  the  specimen  and  secured  by  a  ceramic 
cement.  The  grooves  were  0.56  mm  (0.022  in.)  wide  and  0.5  mm  (0.02  in.) 
deep.  Other  details  of  the  installation  and  procedure  are  given  in 
Reference  1.  The  thermocouple  outputs  were  cross-plotted  to  give  temper¬ 
atures  of  the  mid-chord  at  the  mid-span  at  various  time  increments  after 
immersion  into  the  fluidized  beds.  These  data  are  presented  as  Figure  3 
for  the  five  cases  analyzed.  It  was  assumed  that  there  was  no  temper¬ 
ature  gradient  through  the  thickness  of  the  wedge. 

Another  set  of  thermocouple  data  was  taken  with  five  thermocouples 
mounted  along  the  leading  edge  over  half  the  span.  These  data  revealed 
a  longitudinal  (along  the  span  of  the  wedge)  temperature  gradient  which 
varied  with  the  different  time  increments.  The  maximum  variation  was 
about  16  percent  greater  at  the  ends  of  the  wedge  compared  to  the  mid¬ 
span  and  occurred  after  30  seconds  of  heating.  However,  for  any  one 
time  increment  it  was  found  that  the  ratio  of  the  leading  edge  mid- 
span  temperature  to  that  of  any  other  span  location  was  nominally  the 
same  for  the  five  investigated  cases.  A  least  square's  best  fit  parabola 
was  determined  for  each  time  increment  and  this  is  presented  in  Table  4. 
This  parabolic  temperature  variation  along  the  span  was  assumed  over  the 
complete  chord  of  the  wedge. 

The  temperatures  at  mid-span  were  determined  from  the  appropriate 
plot  in  Figure  3.  For  locations  other  than  mid-span,  the  temperatures 
were  determined  by  using  the  mid-span  temperature  modified  by  the  values 
given  in  Table  4.  Therefore,  the  use  of  Figure  3  and  Table  4  determined 
the  temperature  distribution  at  any  point  of  the  wedge. 
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SECTION  III 

DESCRIPTION  OF  ANALYSIS 

Both  computer  programs  used  three-dimensional  finite-element 
procedures  to  obtain  an  elastic  analysis  of  the  double-edge  wedge 
geometry  specimen.  The  NASTRAN  program  was  used  to  obtain  an  analysis 
only  for  IN  100  alloy  for  the  time  increment  15  seconds  after  immersion 
into  the  heating  bed.  The  IS03DQ  analysis  was  performed  for  the  17 
heating  and  17  cooling  time  increments  (distributed  over  the  3  minute 
immersion  time)  for  each  of  the  five  alloys  as  shown  in  Figure  3. 

The  IS03DQ  program  was  developed  under  contract  by  the  Air  Force  for 
elastic  analysis  -  specifically  aircraft  gas  turbine  blades,  vanes,  and 
disks.  The  NASTRAN  program  was  developed  by  NASA  for  elastic  analysis 
of  generalized  structures.  Documentation  of  the  IS03DQ  program  includes 
a  descriptive  report  (Reference  19)  and  a  user's  manual  (Reference  20). 
Documentation  of  the  NASTRAN  program  includes  a  theoretical  manual 
(Reference  21),  a  programmer's  manual  (Reference  22),  a  user's  manual 
(Reference  23),  and  a  demonstration  problem  manual  (Reference  24).  For 
general  information  on  the  programs,  the  reader  is  referred  to  these 
manuals.  Specific  information  on  how  the  wedge  was  modeled  and  analyzed 
using  these  programs  is  presented  in  the  following  sections. 

1.  IS03DQ  COMPUTER  PROGRAM 

The  model  for  the  double-edge  wedge  was  one-fourth  of  the  structure 
as  shown  in  Figure  1(b).  There  are  reflective  planes  of  symmetry  at  the 
mid-chord  and  mid-span  for  this  structure.  The  nodal  constraints  on  this 
model  (using  the  axis  notation  given  in  Figure  1(b))  are: 

(1)  No  z-displacement  for  nodes  on  the  mid-span  plane  because  of 
reflective  symmetry. 

(2)  No  y-displacement  for  nodes  on  the  mid-chord  plane  because  of 
reflective  symmetry. 

(3)  No  x-displacement  for  the  two  nodes  at  x=0.0  of  the  mid-span 
plan  to  obtain  a  reference  for  displacements. 
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The  coarse  mesh  model  selected  consisted  of  306  nodes  for  the  64 
isoparametric  elements.  A  typical  element  is  shown  in  Figure  1(b). 

The  element  had  mid-point  nodes  along  the  x-direction  but  not  the  y-  and 
z-directions  so  that  each  element  consisted  of  twelve  nodes.  The 
discretization,  including  element  and  nodal  identification,  was  done  using 
a  mesh  generator.  This  pre-processor  (MESH3)  is  part  of  the  IS03DQ 
family  of  programs.  This  program  required  only  the  cross-section  geometry 
of  the  wedge  and  some  mesh  parameters  for  the  geometry  input.  The  maxi¬ 
mum  aspect  ratio  for  the  elements  was  less  than  13. 

Values  for  the  two  temperature  dependent  properties  (modulus  of 
elasticity  and  mean  coefficient  of  thermal  expansion)  were  entered  into 
the  program  as  segments  of  Table  3.  This  table  gives  the  modulus  of 
elasticity  and  mean  coefficient  of  thermal  expansion  for  each  alloy  at 
56°C  (100°F)  temperature  increments.  Six  values  of  modulus  and  thermal 
expansion  for  six  given  temperature  increments  (Table  3)  were  put  into 
the  program.  The  program  selected  the  value  for  the  two  temperature 
dependent  properties  for  each  node  by  using  the  nodal  temperature  to 
linearly  interpolate  within  the  table. 

The  temperature  loading  was  entered  by  means  of  a  temperature  table 
of  13  chord  temperatures  at  four  different  span  locations.  The  program 
assigned  a  temperature  to  each  node  by  weighted  interpolation.  Because 
temperatures  were  assigned  to  nodes  rather  than  elements,  a  straight  line 
gradient  between  adjacent  nodes  was  assumed. 

The  output  selected  from  the  IS03DQ  program  were  the  displacements, 
strains,  and  stresses.  All  of  these  values  were  determined  at  the  node 
points.  This  set  of  data  was  put  on  tape  for  use  by  another  program 
called  PR0UT3.  The  latter  program,  part  of  the  IS03DQ  family,  allows  the 
amount  and  format  of  the  output  to  be  varied  without  requiring  the  complete 
program  to  be  rerun.  Both  the  MESH3  (pre-processor)  and  PR0UT3  programs 
have  plot  capability. 

The  IS03DQ  family  of  programs  were  run  on  the  Wright-Patterson  Air 
Force  Base  CDC  6600  computer.  Plots  (including  Figure  1(b))  were  done 
using  a  Calcomp  on-line  plotter. 
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2.  NASTRAN  COMPUTER  PROGRAM 

The  model  for  the  NASTRAN  analysis  was  similar  to  that  used  for  the 
IS03DQ  program.  One-fourth  of  the  double-edge  wedge  was  used  considering 
mid-chord  and  mid-span  planes  of  symmetry  as  shown  in  Figure  1(c).  The 
nodal  constraints  were  identical  to  those  used  in  the  IS03DQ  program  so 
that  a  valid  comparison  could  be  made.  A  fine  mesh  was  used  in  the 
NASTRAN  analysis  so  that  it  might  be  used  as  the  "baseline1’  for  comparison 
The  model  consisted  of  820  nodes  for  the  354  CHEXA2  (hexahedral)  elements. 
A  typical  element  is  shown  in  Figure  1(c).  The  discretization,  including 
element  and  nodal  identification,  was  done  by  hand  -  no  mesh  generator 
was  used.  The  geometry  was  entered  into  the  computer  program  by  listing 
the  coordinates  of  each  node  point  from  the  origin  as  shown  in  Figure  1(c) 
Elements  were  selected  so  that  the  maximum  aspect  ratio  for  any  element 
was  always  less  than  two. 

The  complete  table  of  temperature  dependent  properties  (modulus  of 
elasticity  and  mean  coefficient  of  thermal  expansion)  for  IN  100  alloy 
was  entered  with  56°C  (100°F)  increments  as  given  in  Table  3.  The  pro¬ 
gram  selected  the  value  for  these  properties  for  each  element  by  using  the 
element  temperature  to  linearly  interpolate  within  this  table.  Since 
temperatures  were  assigned  to  nodes  rather  than  elements,  the  element 
temperature  was  determined  by  the  program  by  averaging  the  eight  nodal 
temperatures.  Since  NASTRAN  does  not  have  the  capability  to  input 
temperatures  by  use  of  equations,  all  temperatures  were  first  hand 
calculated  (using  Figure  3(a)  and  Table  4)  and  then  entered  for  each  node 
point. 

The  output  selected  from  the  NASTRAN  program  were  the  displacements, 
single  point  constraint  forces,  and  stresses.  The  displacements  and 
forces  were  given  at  the  node  points  and  the  stresses  were  given  at  the 
element  centroids.  Stresses  at  the  leading  and  trailing  edges  were 
obtained  by  extrapolation  of  plots  through  the  centroids  of  the  elements. 

This  program  was  run  using  level  16.0  of  NASTRAN  on  a  Univac  1110 
computer.  The  plot  given  in  Figure  1(c)  was  done  on  a  Calcomp  plotter 
using  the  NASTRAN  plot  subroutine. 
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SECTION  IV 

RESULTS  AND  DISCUSSION 

The  results  are  presented  and  discussed  in  three  parts.  First, 
comparison  of  the  IS03DQ  and  NASTRAN  analyses  for  the  check  case  are 
presented.  Second,  results  for  all  five  cases  calculated  by  the  IS03DQ 
program  are  presented  at  the  critical  location.  The  critical  location  was 
taken  as  that  point  on  the  blade  which  had  the  maximum  longitudinal  strain 
range  (algebraic  difference  between  maximum  and  minimum  longitudinal 
strain)  throughout  the  complete  heating  and  cooling  cycle.  This  location 
was  on  the  leading  edge  but  not  at  mid-span  because  of  the  longitudinal 
temperature  gradient.  Lastly,  detailed  computer  plots  for  the  five  cases 
are  presented  at  the  times  of  both  maximum  and  minimum  longitudinal 
strain. 

1.  COMPARISON  OF  IS03DQ  AND  NASTRAN  ANALYSIS 

The  comparison  of  the  analyses  of  the  double-edge  wedge  using  IS03DQ 
with  the  coarse  mesh  model  (Figure  1(b))  and  NASTRAN  with  the  fine  mesh 
model  (Figure  1(c))  is  given  in  Figure  4.  The  comparison  shows  very 
good  agreement.  Both  analyses  were  independently  performed  for  IN  100 
alloy  after  15  seconds  of  fluidized  bed  heating.  This  alloy  and  time 
increment  were  selected  as  being  approximately  the  most  severe  combination 
of  all  those  studied  to  accentuate  any  differences  between  analyses. 

Figure  4(a)  gives  the  normal  x-,  y-,  and  z-displacements  along  the 
leading  and  trailing  edges.  These  results  show  that  the  normal  displace¬ 
ments  as  determined  by  the  two  methods  essentially  coincide. 

Figure  4(b)  gives  the  longitudinal  stress  along  the  mid-chord  at 
one-quarter  span  which  was  the  critical  location  for  this  case.  These 
very  good  comparative  results  show  that  both  the  leading  and  trailing 
edges  are  in  compression.  This  is  due  to  the  manner  of  testing  in  that 
the  specimens  were  stacked  so  that  they  were  heated  and  cooled  from  both 
the  leading  and  trailing  edges.  A  force  balance  of  this  cross-section 
showed  that  equilibrium  requirements  were  satisfied. 
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This  comparison  confirmed  that  the  IS03DQ  program  using  a  coarse 
mesh  model  was  sufficiently  accurate  to  obtain  very  good  quantitative 
results.  It  also  gave  confidence  In  the  use  of  this  specialized  blade 
and  disk  stress  analysis  program. 

2.  CRITICAL  LOCATIONS 

Results  for  the  five  analyzed  cases  at  the  two  critical  locations 
(symmetrical  about  mid-span)  as  a  function  of  time  after  immersion  into 
the  fluidized  beds  are  given  in  Figure  5.  This  figure  shows  the  temper¬ 
ature,  and  longitudinal  strain  and  stress  as  a  function  of  cycle  time 
which  occur  at  both  critical  locations  on  the  leading  edge. 

In  Figure  5,  the  temperature  is  the  nodal  temperature  at  the 
critical  locations  on  the  leading  edge  as  determined  from  the  temperature 
loading  that  was  input  to  the  IS03DQ  program.  The  procedure  used  to 
determine  this  temperature  is  given  in  the  section  IS03DQ  Computer 
Program. 

Both  the  longitudinal  leading  edge  stress  and  strain  show  very  steep 
gradients  for  about  the  first  10  seconds  of  immersion  in  both  the  heating 
or  cooling  beds.  The  results  show  that  the  leading  edge  goes  into 
compression  upon  immersion  into  the  heating  bed.  As  the  specimen  reaches 
a  steady-state  condition,  the  stresses  and  strains  approach  zero.  Upon 
immersion  into  the  cooling  bed,  the  leading  edge  goes  into  tension  fol¬ 
lowed  by  a  gradual  drop-off  to  low  stress  and  strain  by  the  end  of  the 
cooling  cycle. 

The  maximum  longitudinal  strain  range  for  the  cases  analyzed  varied 
from  0.53  to  0.82  percent  (Figure  5).  Mar-M  200  and  Rene  80  demonstrated 
the  highest  strain  range  of  about  0.8  percent.  Mar-M  302  alloy  showed 
the  lowest  strain  range  of  the  five  alloys  analyzed. 

Due  to  symmetry,  the  analysis  showed  three  critical  locations  on 
the  leading  edge  which  were  0.64  cm  (0.25  in.),  1.27  cm  (0.5  in.)  or 
2.54  cm  (1.0  in.)  away  from  mid-span  for  the  five  cases  evaluated. 
Preliminary  experimental  data  In  fluidized  bed  tests  for  some  of  the  cases 
indicate  that  cracks  ire  initiated  in  this  region. 
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The  leading  edge  of  the  double-edge  wedge  is  in  a  uniaxial  state  of 
stress.  On  the  free  surfaces,  the  normal  x-  and  y-stresses  (refer  to 
Figure  1(b)  for  the  axes  convention)  are  zero.  Therefore,  the  effective 
stress  at  the  leading  edge  Is  equal  in  magnitude  to  the  longitudinal 
z-stress.  The  x-  and  y-strains  at  the  leading  edge  equal: 


£ 


X 


-ve_ 


(1) 


where  e  =  strain  in  x-,  y-,  or  z-direction,  and  v  =  Poisson's  ratio. 
By  definition  (Reference  25)  effective  strain  is: 


e2)2  +  (e2  -  e3)2  +  (e3  -  e-j)2 


(2) 


where  1,  2,  and  3  refer  to  the  principal  directions.  Since  the  shear 
strains  are  zero  at  the  leading  edge,  the  normal  strains  equal  the 
principal  strains.  Substituting  Equation  1  In  Equation  2  gives  the 
effective  strain  at  the  leading  edge  as: 


'eff 


=  2(1  +  v) 


(3) 


3.  MINIMUM  AND  MAXIMUM  LONGITUDINAL  STRAIN 

Results  for  the  five  analyzed  alloys  at  the  time  increment  of 
minimum  and  maximum  leading  edge  longitudinal  strain  are  shown  in 
Figures  6  and  7,  respectively.  The  complete  distribution  of  temperature 
and  also  normal,  shear,  and  effective  stresses  and  strains  are  shown  over 
the  complete  mid-chord  plane  of  the  wedge.  The  notation  used  is  con¬ 
ventional  elasticity  notation  with  the  axes  convention  as  given  in 
Figure  1.  The  assumption  of  constant  temperature  through  the  thickness 
of  the  wedge  results  In  zero  y-stress  over  the  mid-chord.  For  this 
reason  the  y-stress  plot  is  not  presented.  The  minimum  (largest  com¬ 
pressive)  longitudinal  strain  always  occurred  during  heating  and  the 
maximum  longitudinal  strain  always  occurred  during  the  cooling  part  of 
the  cycle.  These  plots  were  made  utilizing  the  PR0UT3  program.  These 
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results  show  the  reflective  symmetry  about  mid-span.  These  plots  in 
addition  to  those  in  Figure  5  will  be  used  for  further  evaluation  of 
various  life  prediction  theories  such  as  strain  range  partitioning. 
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SECTION  V 

SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

The  elastic  stress  analyses  for  a  double-edge  wedge  geometry  specimen 
cycled  in  fluidized  beds  were  determined  using  conventional  three-dimensional 
finite-element  techniques.  The  analyses  were  performed  as  a  joint  program 
of  the  Aero  Propulsion  Laboratory  (AFWAL/POTP)  and  the  National  Aero¬ 
nautics  and  Space  Administration  (NASA)  Lewis  Research  Center.  IN  100 
alloy  was  analyzed  using  the  NASTRAN  computer  program  with  a  fine  mesh 
for  only  one  severe  heating  time  increment.  The  Aero  Propulsion  Laboratory 
used  the  IS03DQ  program  with  a  coarse  mesh  model  for  this  combination 
and  all  other  combinations.  Five  alloys  (IN  100,  Mar-M  200,  Mar-M 
302,  NASA  TAZ-8A,  and  Rene  80)  subjected  to  the  same  thermal  cycling 
condition  were  analyzed.  This  condition  was  alternate  3-minute  immersions 
in  fluidized  beds  maintained  at  316°  and  1088°C  (600°  and  1990°F). 

Specific  major  results  are: 

(1)  The  analyses  showed  the  leading  edge  of  the  double-edge  wedge 
goes  into  compression  when  immersed  into  the  heating  bed  followed  by 
tension  when  immersed  into  the  cooling  bed.  Steep  stress  and  strain 
gradients  occurred  during  the  first  10  seconds  of  immersion  in  either  bed. 

For  example,  0.48  percent  strain  was  noted  for  IN  100  alloy  during  the 
initial  5  seconds  immersion  in  the  heating  bed. 

(2)  The  maximum  longitudinal  strain  range  (algebraic  difference 
between  maximum  and  minimum  longitudinal  strain)  for  the  five  alloys 
analyzed  varied  from  0.53  to  0.82  percent. 

(3)  The  two  locations  of  maximum  longitudinal  strain  range  at  the 
leading  edge  of  each  wedge  were  between  0.64  and  2.54  cm  (0.25  and 

1.00  in.)  away  from  mid-span  for  the  five  alloys  analyzed.  Experimental 
test  data  for  the  alloys  that  have  cracked  indicate  that  the  cracks 
initiated  at  these  locations. 

(4)  The  comparison  of  the  analyses  using  a  fine  mesh  model  (354 
elements)  and  NASTRAN  with  a  coarse  mesh  model  (64  elements)  and  IS03DQ 
showed  very  good  agreement  for  the  single  condition  checked. 

(5)  The  results  from  this  Investigation  can  be  used  for  further 
evaluation  of  various  life  prediction  theories. 
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TABLE  1 


ALLOYS  AND  CONDITION  ANALYZED 


Alloy 

Fluidized  bed  cycling  condition  for  all  alloys 

IN  100 

Mar-M  200 

Mar-M  302 

NASA  TAZ-8A 

Rene  80 

Keating  bed  temperature:  1088°  C  (1990°  F) 
Cooling  bed  temperature:  316°  C  (600°  F) 
Immersion  time  in  each  bed:  180  seconds 

TABLE  2 

TEMPERATURE  INDEPENDENT  ALLOY  PROPERTIES 


ALLOY  PROPERTIES 


Alloy 

Poisson's  ratio 

Density 

g/cm'^ 

lb/ in3 

IN  100 

0.2981 

7.750 

0.280 

Mar-M  200 

.  3039 

8.  525 

.308 

Mar-M  302 

.2938 

9.217 

.333 

NASA  TAZ-8A 

.  3166 

8.  636 

.312 

Rene  80 

.3217 

8.166 

.295 
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TABLE  4 

TEMPERATURE  VARIATION  ALONG  SPAN 


2 

[Tx  z  =  Tx  mg  (Az  +  Bz  +  C),  where  Tx  is  the  temperature  at  any  x,  z 

coordinate  (see  fig.  1),  Tx  mg  is  the  temperature  at  the  x  coordinate  at 
midspan,  and  z  is  the  span  coordinate;  all  temperatures  in  °F  (F  = 

9/5  C  +  32).] 

Time  increment, 

Heating  bed 

Cooling  bed 

sec 

A 

B 

C 

A 

B 

C 

0 

-0. 00870 

0.0517 

0.  9205 

-0. 00666 

0.03957 

0. 9427 

3 

.04401 

-.2614 

1.  3891 

-.01775 

.  1055 

.  8447 

G 

.03739 

-.2221 

1.  3290 

-.02384 

.  1416 

.7911 

9 

.  03G88 

-.2191 

1.3372 

-.02548 

.  1514 

.7786 

12 

. 0380G 

- . 22G1 

1.3344 

-.02731 

.  1622 

.7622 

15 

. 03G95 

-.2195 

1.  3300 

-.02889 

.  1716 

.7480 

30 

.02758 

-. 1638 

1.  2504 

-.03047 

.  1810 

.7338 

45 

. 017G9 

-. 1051 

1.  1630 

-.03141 

.  1866 

.7224 

GO 

.01432 

-.08506 

1.  1324 

-.03442 

.  2044 

.6905 

75 

. 0100G 

-.05978 

1.0934 

-.03265 

.  1939 

.7093 

90 

. 00833 

-.04948 

1.0791 

-.02867 

.  1703 

.7440 

105 

.00557 

-.03311 

1.0528 

-.02445 

.  1452 

.7843 

120 

. 00G27 

-.03722 

1.0571 

-.02276 

.  1352 

.7981 

135 

.00440 

-.02614 

1.0415 

-.01876 

.  1142 

.  8323 

150 

.00371 

-.02205 

1.0357 

-.01533 

.09107 

.  8622 

1G5 

. 00297 

-.01762 

1.0285 

-.01278 

.07593 

.  8832 

180 

.00262 

-.01553 

1.0243 

-.01212 

.07198 

.  8876 

15 
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(a)  IN  100  alloy. 

Figure  1  -  Temperature  of  midchord  at  midspan  at  various  times  after  immersion  into  the  fluidized  beds. 


20 


AFWAL-TR-80-201 3 


L 


2000 

1900 

1800 

1700 

1600 

1500 

1400 

1300 

1200 

1100 

1000 

900 

800 

700 

600 


into  heating  bed. 


Figure  3.  -  Continued. 
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Temperature, 


Time  after  immersion 
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Time  after  immersion 
into  heating  bed, 


Position  along  midchord 
le)  Rene'  80  alloy. 
Figure  31  -  Concluded. 
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(a)  Normal  displacements. 

Figure  4.  -  Comparisons  determined  by  using  IS0300  and  NASTRAN  computer  programs  (using  the  models  in  fig.  1)  for  IN  100 
alloy  after  15  seconds  heating  in  the  1088°  C  (1990P  F)  fluidized  bed. 
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Midspan  (plane  of  symmetry) 


IT 


yxy,  10  "mlm  yxz,  10  ’’m/m  y^,  10  ^m/m  ymax,  it 

(a)  IN  100  alloy  after  9  seconds  immersion  in  the  cooling  bed. 


Figure  1.  -  Temperature,  stress,  and  strain  distribution  of  midchord  at  time  of  maximum 
toading  edge  longitudinal  strain.  IF  ■  9/5  C  +  32)11  ksi  *  6. 89x1$  N/nA 


ln.  10 -m/m 

(e)  Rene  80  alloy  after  6  seconds  immersion  in  the  cooling  bed. 
Figure  7.  -  Concluded. 
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